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Abstract 



c/j ■ This paper proposes a novel generalization of group testing, called multi-group testing, 

which relaxes the notion of "testing subset" in group testing to "testing multi-set". 

I> ■ The generalization aims to learn more information of each item to be tested rather 
O ; 

^ ■ than identify only defectives as was done in conventional group testing. This paper 

^ . provides efficient nonadaptive strategies for the multi-group testing problem. The 

^ ■ major tool is a new structure, g-ary additive {w, (i)-disjunct matrix, which is a gener- 
alization of the well-known binary disjunct matrix introduced by Kautz and Singleton 

txj : in 1964. 
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1 Introduction 

A frequently used tool to identify an unknown set of defective (positive) elements 
out of a large collection of elements by group tests is called Group Testing. In the 
classic group testing, a "group" test can be any "subset" of the given collection and 
its outcome is binary under Boolean operations: YES or NO. The former indicates 
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that there is a positive element in this test and the latter implies no positive elements. 
Due to a diversity of its applications, there have been many variants of the classic 
group testing in the literature. Readers are referred to the book [11] and some recent 
papers [1, 4, 5, 6, 7, 8, 16, 19] for further information. 

This paper introduces a novel but natural generalization of Group Testing, called 
multi-group testing. Namely, a test can be applied to any "multi-subset" of the given 
set, where "multi" means every single element is allowed to be taken more than 
once in the same test. Obviously, allowing a number of duplicate copies in a test is 
meaningless under the assumption of Boolean operations. Therefore, in the multi- 
group testing the outcome rule relaxes to arithmetic operations (will be specified 
later). Moreover, the status of each individual item to be tested also relaxes from 
{positive, negative} to a larger set of categories. 

A mathematical model can be described roughly as follows: Let x = {xi,X2, - ■ ■ , Xn) 
e M" be an unknown vector, where Xi denotes the status of ith item. A measurement 
can be apphed to any vector y e {0, 1, • • • }" with an outcome 



The goal is to learn the unknown vector x through measurements in an efficient 
fashion (less measurements and fast decoding). 

An obvious feature under this model is that any measurement which is linearly 
dependent on some other measurements is useless. The reason is that its outcome 
can then be simply derived from a linear combination of the outcomes of the others. 
The other feature is that, without any further information of the unknown vector x, 
n measurements are necessary in the worst case to learn the unknown vector (n is 
clearly sufficient) . The reason is that after k measurements x can be any vector in the 
n — k dimensional subspace whose outcome is consistent with the k measurements. 
If A; < n, such vectors are not unique and thus cannot be determined exactly. 

To make the multi-group testing model more interesting and challenging, we shall 
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assume that 

1. the unknown vector x is d-sparse, that is, x contains at most d nonzero entries 
where d is a constant with d n; 

2. the entries of x all belong to a certain set D C M which is a priori knowledge; 

3. the number of copies from every single item in a measurement is restricted to 
{0, 1, • • • ,q — 1}, i.e., y e Z^, where the integer q is prescribed. 

Note that the outcome Vy — Yli=i Vi^i possibly exceed q. 

This paper focuses on nonadaptive strategies where measurements are performed 
simultaneously and therefore all measurements must be settled in advance. A non- 
adaptive strategy that uses t measurements can be represented by a ^ x n matrix 
A — [ttij] with columns as items and rows as measurements, and the value at de- 
notes the number of copies of item j in measurement i. The nonadaptive multi-group 
testing problem can then be converted into the problem: Construct a matrix A = [aij] 
with ttij e Zq so that the unknown (i-sparse vector x e can be determined exactly 
and efficiently through the outcome vector v = Ax. 

Motivation and related work 

This generalization naturally arises in several situations. Clearly, the conventional 
quantitative group testing in [10] is a special case of the multi-group testing with 
D = {0, 1} and g = 2. 

In some applications, such as blood testing, what patients or doctors want to 
know might be not only a yes-or-no answer but also a more precise index, a standard 
by which the level of some illness can be judged. The major purpose of relaxing x 
from the usual set {0, 1} to a prescribed set D (can be very large) is that we aim 
to learn more information that each item carries rather than just determine which 
items are positive. Of particular note is that D is allowed to contain not only positive 
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elements but also negative elements, whose presence is in a sense to cancel the effect 
of positive elements; such a phenomenon was noticed in real applications: molecular 
biology [14], blood testing [17] and drug discovery [18], and discussed in the hterature 
of group testing as the inhibitor model [3]. 

The multi-group testing problem is also closely related to the compressed sensing 
theory, a popular and intensively studied field in the past decade (see [2, 9] and ref- 
erences therein). Compressed sensing is a signal processing technique for recovering 
a signal by finding solutions to underdetermined linear systems (more unknown vari- 
ables than equations). Typically, the signal to be recovered is assumed to be sparse, 
which coincides with the multi-group testing problem. However, in multi-group test- 
ing the measurements are restricted to vectors in while in compressed sensing they 
are usually vectors in the real space R". A consequence of the fundamental difference 
is that the developed theory in compressed sensing cannot be applied immediately to 
the multi- group testing problem because methods work for M need not work for Z^. 

The "multiset" relaxation has an advantage of performance. Using less measure- 
ments in the multiset model than in the set model is to be expected as set is a special 
case of multiset. To the best of our knowledge, there are a few studies about non- 
binary structures relevant to group testing. However, nonbinary structures are used 
only to construct binary designs for group testing purpose in almost all the studies 
except Emad and Milenkovic [13] , where a semi-quantitative model is considered. 

Our contribution 

We give nonadaptive strategies for the multi-group testing problem with general D. 
We note that the one-sided case, i.e., D is nonnegative or nonpositive, is a relatively 
simple case to handle. The reason is that in this case a zero outcome simply implies 
that all the items appearing in the measurement are zero, in contrast to the general 
case, a zero outcome can be produced by a combination of some positive elements 
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and negative elements. 

Our decoding algorithms are quite efficient. For the one-sided case our strategy 
has a decoding algorithm in 0{\D\tn) time and for the general case it is 0(|D|tn'^+^). 
By contrast, the brute-force method needs (^\D\Hn = 0{\D\Hn'^^^) time in the 
worst case. 

The rest of the paper is organized as follows. Section 2 introduces our major tools 
and studies the one-sided case. Section 3 deals with the general case. Finally, Section 
4 provides some explicit solutions to the multi-group testing problem. 

2 The one-sided case 

This section starts with our major tools. 

Definition 1 A q-ary matrix M — [m^] of size txn is called additive {w, d)-disjunct 
if for any s G {(si, S2, - ■ ■ , Sn) G {0, l}" : ^ d,} and for each k G [n] such that 

Sk — 0, there exists an i & [t] such that 

n 

Remark: By definition, the well-studied d-disjunct matrix [12, 15] (also known as 
cover-free families or superimposed codes) is merely a binary {q = 2) additive {w, d)- 
disjunct matrix. In particular, when w > q, any g-ary additive (w, (i)-disjunct matrix 
can be converted simply to a binary d-disjunct matrix by replacing every non-zero 
entry with 1. 

Theorem 2.1 Let A be a matrix in R*^". Given a fixed (unknown) vector x G M", 
let V = ylx. Let x' = ax -|- b for some ^ a E M. and b G M" (known), and let 
v' = A'x'. Then the problem of learning x from v is equivalent to the problem of 
learning x' from v'. 
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Proof. The proof follows immediately from the linear mapping from v' to v, that 
is, v' = Ax.' = A(ax + b) = aAx + Ah ^ av + Ah. ■ 

To present our algorithms, we first introduce some notations. Throughout this 
paper, let x e be an unknown d-sparse vector, A — [aij] be the matrix corre- 
sponding to the measurements and v = Ax be the outcome vector. For any vector 
y = iVu ■ ■ ■ ,yn): define 

0(y)^|{y,:y,^O}|. 
Given any vector y = (yi, • • • ,yt), for each j G [n] define 

tj{y) = \{i: ttij >yi}\. 

For any vector y = {yi, • • • , yn) and 5 e R, define 8^'='^ = (•^i'^'^j ^2'^^^ ' ' ' > ^n'^^) 
where 

1 otherwise. 

For convenience, we shall use Sj to denote if no confusion occurs without specified 
superscripts. Likewise, we define s^'-'^ and s^'-'^ by replacing yj — S with yj > S and 
yj < respectively. For any two vectors x and y of the same dimension, denote by 
X ^ y if > i/i for all i and by x 5:^ y otherwise. 

The rest of this section focuses on the one-sided case. By symmetry, we may and 
shall assume that D = {0, Ci, ■ • ■ , c^} where < Ci < C2 < • • • < c^. 

Lemma 2.2 Let x e {0,ci, ■ ■ ■ be an unknown sparse vector with (/){x) < d, 

where < ci < C2 < • • • < c^. Let A he a q-ary additive {w, d)-disjunct matrix of size 
t X n with w > — and v = Ax be the outcome vector. Then s^'-'^^ can be identified 

— Cl •' 

from V. 

Proof. By Theorem 2.1, this problem can be reduced to the problem of learning the 
unknown sparse vector x' = (x'l, x^^ ■ ■ , x'^) G {0, 1, ^, • • • , from the outcome 
vector v' = —v. 

Cl 



Cm/Cl / 

For each k such that x'j. 7^ 0, for all i e [t] we have aik < ^ c ^ aij | = u; . 

Accordingly, tk{v') — whenever x'j^ ^ 0. 

Since 0(x) < c?, '-^ < d. For each k such that x'^ = 0, by definition of the 

g-ary additive (w, (i)-disjunct matrix, there exists i e [t] such that 



j=i j=i 

This implies that tkiv') > if = 0. By the above discussion, we can identify the 
vector s^''-^ through the counting function ife(v'). As a consequence, s^'-'^^ can be 
identified too. ■ 

We now analyze the time complexity for the decoding algorithm corresponding 
to Lemma 2.2. For each j G [n], it takes t operations for computing the value tj{-). 
Therefore, the decoding complexity is 0{tn). 

Corollary 2.3 Let x e {cq = 0,ci,--- ,0^}" be an unknown sparse vector with 
(/»(x) < d, where < Ci < C2 < • • • < c^. Let w — maxi<fc<^{ "^Zll'l } • Then 
any q-ary additive {w, d)- disjunct matrix of size t x n can be used to identify x with 
0{\D\tn) decoding complexity. 

Proof. The corollary follows by applying Lemma 2.2 repeatedly (m times). The 
precise process is as follows. Let x^ = ^x and = (here we shall use x' and 
v' to denote the updated vectors in the ith round). Since w > by Lemma 2.2 we 
know that s''^'-^ can be identified from successfully. 
For k — 2, - ■ ■ ,m — 1, define recursively that 



_ ^ (v- - ^^"•>-^) and x' = ^ (x'-' - ^■'•-■i') . 



Cfc— Cfc-2 

Cfc-1— Cfc_2 ~ Cfc_i— Cfc_2 



It is easily verified that v'^ = Axl" for all k. Note that x*^ e {0, 1, ^i+iz£h^ ... ^"-'^"-n " 
is a (i-sparse vector and w > '^^Z^l'l for each A; = 1, ■ ■ ■ ,m. Applying Lemma 2.2 



repeatedly, '-^ can be identified for all /c = 1, ■ ■ ■ , m. Consequently, the unknown 
vector X can be identified as x = Yl^i (^i — Ci-i)s^''-^ — X^^i QS^'"'^*. This com- 
pletes the proof. ■ 

3 The general case 

In this section, we turn our attention to the general case that elements in D are 
neither all nonnegative nor all nonpositive. Let D — {zm2,-'' > -^i, 0, ci, • • • ,c^i}, 
where z^^ <---<-2;i<0<ci<---< Cmi- For any vector v, r e {0, 1}"^ and /i e R, 
define 

v'*'" = V + hAr. 

For any vector v, o? e N and /i e M, define 

t*(v, h, d) = min tjiv'^'') for j — 1- • ■ ,n. 

Theorem 3.1 Let x e {zm2, • • • > -^i, 0, Ci, • • • , c^i}", Zm2 < • • • < -Zi < < Ci < 

• • • < Cmi Ojud (/)(x) < d, he an unknown sparse vector. Let w > £!i!i_f!ii2 Then any 
q-ary additive {w, 2d) -disjunct matrix of size t xn can be used to identify s'''-'^^ 

Proof. By Theorem 2.1, this problem is equivalent to the problem of learning x e 
{^, • • • , a, 0, 1, |, • • • , ^}-. Let X = (xi, • • • , A be a g-ary additive {w, 2d)- 
disjunct matrix, w > Si^iZfuis.^ of size t x n and v = Ax be the outcome vector 



^"12 ... 21 Q ][ £2 . . . £ 

CI ' ' CI ' ' ' CI ' ' CI 

Cl 

corresponding to x. It suffices to identify s'''-^ 



Consider any vector r(x) = (ri,r2, • • • ,r„) e {0, 1}" with 0(r(x)) = d such that 
r(x) y s^'^°. For each k with > 1, for all i e [t] we have 

c, \ ''m2 \ / . ^r(x) -^1712 4 / \\ 

^ — Vi > Oij-rj (smce v ^=1 ^ — v ^^\^)) 

n 

^ ^ + ^ij ~ ^^^^^ (since r(x) ^ s'''^^) 

{i:rr,>l} {j-Xj<0} j=l 

{i:a:,>l} 



By definition, tk{v ~'^^^') = and therefore t^(v, — d) = if Xk > 1. 



Cl 



Consider tfie case Xk < 1. For any arbitrary vector r e {0, 1}" witfi 0(r) < d, 
consider the unknown vector s = (si, S2, • • ■ , s„) where 



1 if > 1 or = 1, 
otherwise. 



As 0(r) < d and 0(x) < d, we have 0(s) < 2d. Since A is a g-ary additive (w, 2d)- 
disjunct matrix where w > ^""i ^ ^j^ere exists i G [t] such that 



Oife > — OySj ( by definition of disjunctness) 



n n 

— aijSj — aij 

n n 



— ^ij-^'^j ~ '^ij^j ( since s >: x) 

n 



n 

Cl ^ ' ' 

n 



> Vi — aijfj ( since s ^ r(x) and < 0) 



This imphes ife(v~^'') > for any r e {0, 1}" with 0(r) < d. Hence, i^(v, d) > 
if Xfc < 1. 

Therefore, by the above discussion, one can determine whether Xk > 1 through 

the counting function tl{v, — d). ■ 

We now analyze the time complexity for the decoding algorithm corresponding to 
Lemma 3.1. For each j e [n], it takes t(^) operations for computing the value tj{-). 
Therefore, the decoding complexity is 0{tn'^'^^). 

Combining Corollary 2.3 and Theorem 3.1, we have the following result. 

Theorem 3.2 Let x e {2:^125 ' ' ' ,zi,0,ci, - ■ ■ , c^i}", < ■ ■ ■ < zi < zq = = 
Co < Cl < • • • < c^^ and 0(x) < d, be an unknown sparse vector. Let w > 



max 



Cm I Ci Z< 



m2 



0<i<mi — ill J< — '^'"2 + 1 ]^ < ^ < ^2 I I . Then any 

Ci+l - Ci } ^ { + ) ) 

q-ary additive {w, 2d) -disjunct matrix of size t x n can be used to identify x with 

0{\D\tn''-~^^) decoding complexity. 

Proof. By Theorem 3.1, s^'-'^i can be identified since w > ^iiH—fziis. = x — 

^^gx,>ci g^j^j for i 2, • • • , mi — 1 define recursively x* = x'~^ — (q — Ci^i)s^'-'^' . Note 
that x' e {zjn2, • • • , -Zi, 0, Cj+i — Cj, • • • , Cjni — Cj}" and (/>(x*) < for i = 1, • • • , mi — 1. 
With w > ^""^ ^"'^ , applying Theorem 3.1 repeatedly, one can successfully identify 
sx%>c,+i-ci j-Qj. equivalently s^'-'^'+i) for i = 1, • • • , mi — 1. As a result, one can identify 
gx,=ci _ gx,>ci _ gx,>ci+i for alH = 1, • • • , mi (for well-defineness let c^j+i = oo). 

( mi \ 
^ Qs'^'^'^* - X j e {0,-2;i,--- ,-Zm2T] hence 0(x') < d. As 

w > max |^!!^2_lL^!zi\ by Corollary 2.3, x' can be identified. This completes the 

l<i<m2 —Zi + Zj-i 

proof. I 

Note that the bound on w in Theorem 3.2 is not necessary the best choice. One 
might obtain a better bound by first applying the transformation method introduced 
in Section 2. The following demonstrates such an example. 

Example 1 Letx^ e {-2,0,1,4}" with < d and let = -x^ e {-4,-1,0,2}'^ 

By Theorem 2.1, we know that the problem of learning x-*^ is equivalent to the prob- 
lem of learning x^. However, to apply Theorem 3.2 to x^, we have to require vu'^ > 
max{^^3^, — 4 which is smaller than 6 = max{ ^~^^^ , , < w^, 

required by simply applying Theorem 3.2 to k'-. ■ 

4 Constructions for multi-group testing 

This section proposes three constructions for the multi-group testing problem. The 
first one is to construct a g-ary additive {w, d)-disjunct matrix from a conventional 
binary disjunct matrix by deleting some rows. 
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Theorem 4.1 Let A be a q-ary additive {w,d)- disjunct matrix of size t x n and 
w < If A has d+2 rows Ri, • ■ ■ , Rd+2 pairwise disjoint with entries in {0, 1}, then 
A' obtained from A by deleting Rd+2 o-nd replacing R^ with R'^ — {wd-\- l)Rk + Rd+2 
for k = 1, ■ ■ ■ , d + 1 is a q-ary additive {w, d)-disjunct matrix of size {t — 1) x n. 

Proof. For convenience, represent A using row indices [t] and column indices [n] 
and without loss of generality assume Ri,--- ,Rd+2 be the first d + 2 rows, i.e., 
indexed from 1 to d + 2. Let A' be the obtained matrix (using the same indices 
with A). Obviously, every entry in A' is at most q. Consider fixed j e [n] and 
Ji, • ■ ■ ) Jd £ W \ {j}- Since A is q-ary additive {w, (i)-disjunct, there exists an i e [t] 
such that Uij > w Ylt=i '^Vk- There are only three cases as follows. 
U i ^ {!,■■■ ,d + 2}, then row i is in A' and a'^j = aij > w ^^^^ a^j^ = w Yfk^i o^^, 
as desired. 

If i e {1, • • • , (i + 1}, then row i is in ^4', = 1 and aij^, = for A; = 1 • • • , 

Therefore, a[j = {wd + 1) > w Yl'k=i ^ ^ Yli'k=i '^ijk '^^lere the last inequality holds 
for a^^ < 1. 

U i — d + 2, then row i is not in A'. We need to find another row i' in A' with the 
desired property. In this case, ad+2,j = 1 and ad+2,j,, = for A; = 1 • • • , (i. Since 

• • • , Rd+i are pairwise disjoint, at least one of them, say i', has all entries at 
the columns j^s for A; = 1, • • • , <i, i.e., a^ij^^ = for A; = 1, • • • , d. Since ad+2,j — 1 and 
therefore Oj/j = 0, we have a^,^- = 1. Thus, a[,j > w Yl'k=i '^'i'jk ~ ^■ 

Since j,ji, ■ • ■ ,jd are chosen arbitrarily, the proof is complete. ■ 

A binary matrix A is called transversal if its rows can be divided into disjoint 
families such that rows in each family are disjoint. We say a family of size b if it has 
b rows. Denote fb{A) as the number of families of size at least b in the matrix A. 
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12x16 

It is easily verified that A is 2- disjunct and transversal as it can be divided into 3 
disjoint families. Further, /4(A) = 3. ■ 

Corollary 4.2 Let A be a transversal d-disjunct matrix of size t x n. There exists a 
q-ary additive {w, d) -disjunct matrix of size t' xn withq > wd-\-2 andt' — t — fd+2{A). 

Proof. The Corollary follows immediately from Theorem 4.1. ■ 
Example 3 Let 
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obtained from A by using operations as in Theorem 4-1- Then A' is a Q-ary (2,2)- 
disjunct matrix of size 9 x 16. I 



Next, we focus on the particular case that D is one-sided. The following two con- 
structions are based on a special operation of matrices, referred to Kronecker Product. 
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Of particular note is that the obtained structures also admit efficient decoding algo- 
rithms. 



Definition 2 If A and B are matrices of size n x m and s x t respectively, then the 
Kronecker Product A® B of the two matrices A and B is the ns x mt matrix whose 
entries {A ® B)ij^ki = Ai^Bji with row indices listed as 11, - ■ ■ , Is, 21, • • • , nl, ■ ■ ■ ,ns 
and column indices as 11, • • • , It, 21, • • • , ml, • • • , mt. 



The Kronecker Product of two matrices A and B can also be viewed as 



A®B 



( a\\B ayiB 

021 -B 022 -B 



yiniB a„2-B ••• anmB j 

Theorem 4.3 Let A he a matrix of size t x n that can successfully identify any 
unknown vector y e {0, Ci, C2, • • • , c^}" with 0(y) < d, where < Ci < • • • < c^. Let 
B he a hinary d-disjunct matrix of size t' x n' . Then B ® A can successfully identify 
any unknown vector x e {0, ci, C2, ■ ■ ■ , Cm}^^' with 0(x) < d. 

Proof. Partition the unknown vector x equally into n' pieces xi, X2, • • ■ , x„/ where 
Xj e {0, ci, C2, • • • , Cm}"' for all j = 1, • • • , n'. Then 



[B ® A)x 



/ buA buA 
b2iA 622^ 



bi„'A\ 
h^n'A 



\ht'iA bt'2A ■■■ bt'n'Aj 



/xi\ 

X2 



V2 



where the outcome vector v = 



/vi\ 

V2 



with each Vj a vector of length t. We now want 



to show that the unknown vector x can be identified from the outcome vector v. 

We first show that one can determine whether x^ = 0, where denotes a vector 
whose entries are all zero, for all j through the d-disjunctness property. Obviously, 
there are at most d x^'s such that Xj 7^ as 0(x) < d. Let J be such a set consisting 
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of indices of these x/s with Xj ^ 0. For a fixed k with x^ = 0, by the definition of 
d-disjunctness, there exists an i such that 6^^ = 1 and % = for all j £ J ■ It follows 





that V£ = i^t\A be2A 



X2 



. In contrast, for each j e J, 



for each z with bij = 1 we have Vj = {bnA bi2A 



bin' A) 



/xA 







since 



fcjj Axj and entries in bij/Axj/ are all nonnegative. As a result, one can determine 
whether Xj is or not by checking Vi — for some i with bij — 1. 

For each j with x^- 7^ 0, by the definition of (i-disjunctness, there exists an i such 
that bij — 1 and bij> — for all j' e J \ {j}. It follows that 

Vi = (6a-4 bi2A ■ ■ ■ bin' A) 



X2 



Ax,-. 



\Xn'/ 



Consequently, Xj can be identified from the outcome Vj since A can successfully 
identify any unknown vector y e {0, ci, C2, • • • , c^}" with 0(y) < d, where < ci < 
• • • < Cm- By the above discussion, one can successfully identify x by using the matrix 
5(g) A I 

Following the idea in [13] of concatenating several matrices, we obtain the fol- 
lowing result: Let A be a binary d-disjunct matrix of size t x n. Let u — Cmd and 
n' — [logy {1 + {u — l){q — 1))J . Construct a g-ary matrix C of size txnn' by concate- 

nating n' matrices: C = (Ai A2 ■ ■ ■ An') , where Aj = j J A for 1 < j < n'; 

\ i=o J 

equivalently C — B ^ A where 



B=(l 1 + u 1 + u + u^ ■■■ ^) . 



Theorem 4.4 Let C be the q-ary matrix as defined above. Then C can successfully 
identify any unknown vector x e {0, Ci, • • • , Cm}""' with 0(x) < d and < Ci < • • • < 
Cm where all Ci 's are positive integers. 
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Proof. Let v = Cx be the outcome vector. Consider x = (xi X2 • • 
Xj's are vectors of length n. We prove this theorem by showing x^j^, 
identified successfully one by one. 
Observe that 



x„') where 
, Xi can be 



- 1 , \ - 1 , 

2^ - — r^^j I + -Axn'. 



U 



u 



(1) 



n'-l 



1 



U 



Since < d and is the largest value, each entry of the term j —Ak 

. . . - - , u — 1 

m (Ij IS at most 



M — 1 y — 1 

This imphes x„/ can be identified from the outcome 



It follows that 
u-1 



Ax„/. 



- 1 

as A is a (i-disjunct 



M"' - 1 

matrix, which guarantees the identification of any d-sparse vector with entries all 
nonnegative. 

Assume now that x„/, x„/_i, • • • ,Xjt+i have been identified. Define = v — 

u—1 

J2j=k+i^j-^r analogous argument shows that 



u'' — 1 



d-disjunctness property again, one can identify x^ from the outcome 



Axfc. By the 

u-1 



— 1 



Vfc 



Repeating this process, x^/, 
completes the proof. 



,xi can be identified successfully one by one. This 
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